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Abstract 

In a recent work (Borras et al, Phys. Rev. A 79, 022108 (2009)), we have determined, for various 
decoherence channels, four-qubit initial states exhibiting the most robust possible entanglement. 
Here we explore some geometrical features of the trajectories in state space generated by the 
decoherence process, connecting the initially robust pure state with the completely decohered 
mixed state obtained at the end of the evolution. We characterize these trajectories by recourse to 
the distance between the concomitant time dependent mixed state and different reference states. 



I. INTRODUCTION 



Entanglement and decoherence are two closely related phenomena that lie at the core of 
quantum physics l|-|3|. Entanglement is nowadays regarded as the most distinctive feature 
of quantum mechanics, and in recent years it has been the subject of intense and increasing 
research efforts. The phenomenon of decoherence consists, basically, of a set of effects 
arising from the interaction (and concomitant entanglement development) between quantum 
systems and their environments 2|, |3|. Every physical system is immersed in an environment 
and interacts with it in some way. The effect due to this interaction is one of the main 
obstacles to the practical implementation of quantum technologies based on the controlled 
manipulation of entangled states such as quantum computing 2|. The decoherence process 
leads the system from a pure state to a (usually less entangled) mixed state. This decay 
of entanglement has recently attracted the interest of many researchers j4-15|. It has also 
been shown that in some cases entanglement can vanish in finite times. This phenomenon 
is known as entanglement sudden death and it has been the focus of numerous contributions 
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-|2l|. This abrupt disappearance o 



entanglement between the reservoirs 



entanglement is closely related to the sudden birth of 



22|. 



Recently, several dynamical properties of entanglement, like asymptotic birth of entangle- 
ment and entanglement sudden death, have been discussed from a geometrical point of view 
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24j . These geometrical interpretations allow for the explanation of the necessary and 



sufficient conditions for the last phenomenon to take place [23|. Various examples according 
to different possibilities for the geometrical details of the set of time asymptotic states are 



provided in 



24|. 



We have recently studied the decay of entanglement under different paradigmatic noisy 
channels, identifying the initial states exhibiting the most robust entanglement j4|. 

The aim of the present contribution is to explore characteristic traits of those state space 
trajectories associated to robust states, with emphasis on the study, from a geometrical point 
of view, of optimal-ones. This will hopefully shed some light on the existence and behaviour 
of robust states. The paper is organized as follows. In Section II we briefly review the local 
decoherence models for multi-qubit systems and the entanglement and distance measures 
that will be used in the present work. In Section III we investigate the evolution trajectories 
in state space of highly entangled four-qubit states. Finally, some conclusions are drawn in 
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Section IV. 



II. PRELIMINARIES 

The systems under consideration in the present study consist of an array of independent 
qubits initially entangled due to a previous, arbitrary interacting process. Each qubit in the 
composite system is coupled to its own environment. In such a local-environment formulation 
there is no communication among qubits and the entanglement between the subsystems 
cannot increase because of the locality of the involved operations. We assume that all 
qubits are affected by an identical decoherence process. The dynamics of any of these 
qubits is governed by a master equation from which one can obtain a completely positive 
trace-preserving map e which describes the corresponding evolution: pi{t) = epi{0). In 
the Born-Markov approximation these maps (or channels) can be described using its Kraus 
representation 

M 

e^P^{0) = J2Ej^p^{0) Ej,, (l) 

i=i 

where Ej j = 1, . . . , M are the so called Kraus operators, M being the number of operators 



needed to completely characterize the channel 



25| . Using the Kraus operators formalism it 



is possible to describe the evolution of the entire N-qubit system, 

p(t) = e p(0) = ^Eii® EjN p(0) [Eii®...®Ej n]'^ . (2) 



A. Decoherence models 

We concentrate our efforts in the study of the following family of decoherence channels: 
the bit flip (BF), phase flip (PF), and bit-phase flip (BPF). These channels represent all the 
possible errors in quantum computation, the usual bit flip 0^1, the phase flip, and the 
combination of both, bit-phase flip. The corresponding pair of Kraus operators {Eq, Ei} for 
each channel is given by: 

Eo = Vl-p/2/, E{ = ^J^Gi- (3) 
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where / is the identity matrix, ai are the Pauh matrices and i = x give us the bit flip, i 



the phase flip, and i = y the bit-phase flip. Following Salles et al. the factor 2 in Eq. 
(E]) guarantees that at p = 1 the ignorance about the occurrence of an error is maximal, and 
as a consequence, the information about the state is minimum. The results obtained with 
the phase damping channel are the same than those of the phase flip channel, actually they 
can be shown to represent the same process, under a proper change of variables 2|. 

We also consider the depolarizing channel (D) which can be viewed as a process in which 
the initial state is mixed with a source of white noise with probability p. For a rf- dimensional 
quantum system, it can be expressed as 

ep = ll+{l-p)p (4) 
The Kraus operators for this process, including all Pauli matrices are 



= ^/l-p'I■, = \l^a, (5) 



where = -j. Under this process the state turns separable after a flnite time, being this 
channel the only one considered in this work exhibiting the phenomenon of entanglement 
sudden death. We remember here that due to the high symmetry of the depolarizing channel 
the evolution of the entanglement depends only on the amount of entanglement of the initial 
state. Then, according to this process, we obtain equivalent evolutions of the entanglement 
for initial states which are equivalent under local unitary operation jj] . Some results for this 
channel will be commented at the end of Sec. III. 

B. Multipartite entanglement quantification 

One of the most popular measures proposed to quantify multipartite entanglement is 
based on the use of a bipartite measure, which is averaged over all possible bipartitions of 
the system. It is mathematically expressed in the fashion 



[N/2] 



E = — V E^""\ (6) 

I- ' J ■m=l 



Here, E{i) stands for the entanglement associated with one single bipartition of the A^-qubit 
system. The quantity gives the average entanglement between subsets of m qubits 
and the remaining N — m qubits that constitute the system. The average is performed 
over the ensuing N^^^^^ nonequivalent bipartitions. If one uses the linear entropy Sl of 
the reduced density matrix of the smaller bipartition to compute E{i), Ef^ turns out to 
be the well known Meyer- Wallach multipartite entanglement measure 26|]. This measure 
was later generalized by Scott to the case where all possible bipartitions of the system were 
considered 27|]. The Meyer- Wallach multipartite entanglement measure has recently been 



related to the regularized quantum Fisher information measure which gives the estimation 
of the strength of low-noise locally depolarizing channels 28 



We will use the negativity as our bipartite measure of entanglement because we are dealing 
with mixed states. The negativity is proportional to the sum of the negative eigenvalues ai 
of the partial transpose matrix associated with a given bipartition. The properly normalized 
negativity reads 



Neg 



1^ 



ai 



(8) 



Distance measures 



Distance measures between quantum states constitute important tools in quantum in- 
formation theory 29|-|33| . In the present work, in order to characterize the trajectories of 
decohered states, we compute the distances between the state of interest and several refer- 
ence states, such as the initial, final, and maximally mixed ones. We measure the distance 
between two different quantum mixed states by recourse to the quantum Jensen Shannon 
divergence (QJSD), which can be defined in terms of the relative entropy as 29 1 



Sip, 



p + cr 



S 



cr. 



p + cr 



(9) 



2 / V 2 

that can be recast in terms of the von Neumann entropy H]\f{p) = —Tr(p log p) in the fashion 



djsip, cr) = 



N 



(10) 



The metric character of the square root of the QJSD has been ascertained recently for 
pure states, and strong numerical evidences have also been found for the mixed states case 



30l . l3l| . To avoid an exclusive dependence on the above distance measure we also use the 
Hilbert-Schmidt distance [32 1. 



dusip^cr) = Wp-crfffs = Tr[{p-af]. (11) 

We will be able to appreciate the fact that the results obtained with both distance measures 
are qualitatively the same. 

D. Robust Maximally Entangled Four-Qubit States 

In a previous work we have performed an iterative numerical search based on the simulat- 
ing annealing algorithm and determined the initial states that evolve (via the just described 



noisy channels) to mixed states with the larger possible amount of entanglement j4|. The 
proof by Higuchi and Sudbery that a four-qubit pure state exhibiting the theoretically max- 
imum amount of entanglement (that is, having all its marginal density matrices maxi mally 
mixed) does not exist constituted a landmark in the study of multiqubit entanglement |34l |. 



In Ref. 



34j . a promising candidate for achieving the maximally entangled state status was 



also proposed, namely. 



l^tot) = \HS) = 11100) + 10011)+ (12) 



(^|1001) + |0110)j +a;2(^|1010) + |0101)j 

with CO = — _| + Such conjectures have later received support from several numerical 
studies |35l-l37|. In our previous work, the \HS) was found to be a robust state. The 
concept of robustness used in this work can be easily explained from Fig. [1] The \HS) state 
is considered to be the most robust state because, for any value of p, it leads to decohered 
states exhibiting more entanglement than those associated with other initial states. In 
this graph we plot the entanglement-evolution of the \HS) state under the action of the 
BF channel and compare it to the entanglement decay of the well-known entangled states 
\GHZ) and \W) 
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FIG. 1: Entanglement evolution for several 4 qubits states under the action of the BF channel. All 
shown quantities are dimensionless. 



\GHZ) = -L(|0000) + 11111)), (13) 
\w) = ^(loooi) + lOOlO) + lOlOO) + |1000)). (14) 

We also found that, for six qubit systems the robust state I'^rob) turns out to be precisely the 
maximally entangled state encountered by some of the authors of a previous work jssl. For 
five qubit systems, the robust state I'^rob) ^^at we found is not as good as the one detected in 
the 4 or 6 qubits instance. For BF and BPF channels its entanglement becomes lower than 



that of other states for large p values {4]. The entanglement decays of I'^rob) ^^"^ l^rob) 
quite similar and their entanglements are always larger than that of any other state tested 
in our samplings. 

III. OPTIMAL TRAJECTORIES 

In order to compare the behaviour of the entanglement decay for the different channels 
we compute it in terms of the degree of mixedness, given by the linear entropy of the density 
matrix of the (evolved) system 
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FIG. 2: Entanglement evolution as a function of the linear entropy of 4 qubits representative 
states under phase flip, bit flip, and bit phase flip decoherence models. All depicted quantities are 
dimensionless. 



Sl{p) = j^il-Trlp']), (15) 

where = 2" and n is the number of qubits. Here, we display only the results corresponding 
to four qubit systems, although similar results are obtained for systems with different number 
of qubits. 

We study the decay of entanglement for different initial states undergoing several decoher- 
ence processes. We note that, for those channels for which a given state is robust, the decay 
of entanglement in terms of the degree of mixedness is the same. In Fig. [2] we plot the decay 
of entanglement for different four qubits initial states for BF, PF, and BPF decoherence 
channels. The entanglement evolution for the \HS) state coincides for the three considered 
channels. Remember that the \HS) state was found to be robust under the action of the 
three channels studied in this work. In contraposition, the entanglement evolution of the 
\GHZ) and \W) states under the PF channel is not equivalent to that under BF and BPF 
processes because these states are not robust for the former development. 

We also compute the distances between the states (at any given time) and some reference 
states using the Q JSD. We considered as reference states the initial, final, and the maximally 
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FIG. 3: QJSD between the decohered state (mitiahy l^'^gf,)) and reference states: maximally 
mixed state (solid line), initial robust state (dashed line), and final separable state (dotted line). 
All plotted quantities are dimensionless. 

mixed state (MM). In Fig. Owe plot the resultant distances between i) the mixed state 
obtained at each time step of the evolution of the inial robust pure state of 4 qubits and ii) 
the reference states for the previously mentioned decoherence processes. Similar coincident 
curves (not shown) are obtained for \GHZ) and \W) states if we consider only the BF 
and BPF channels. According to both graphs the state with robust entanglement (w.r.t. 
several decoherence channels) apparently evolves in the same manner under the action of 
such channels. It is important to note that these trajectories are not actually the same, they 
are just equivalent. 

Finally, we compute the distances between the three final states. These resulting states 
generated by the PF, BF, and BPF channels when acting upon the initial state |$rofe) turn 
out to be equidistantly distributed, i.e., the distance between any pair of them is always 
the same. The distance from any of them to the maximally mixed state or to the initial 
state is also always the same. The overall picture is displayed in Fig. |H The initial state 
is represented by the black square, and the three final states are represented by the black 
spots. These final states are placed at the border of the set of separable states, represented 
by the grey sphere. The star placed in the middle of the sphere denotes the maximally mixed 
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FIG. 4: Equivalent trajectories in Hilbert space, the initial robust 4 qubits state is represented by 
a black square and the final states corresponding to different decoherence processes are represented 
by black spots at the border of the set of separable states (grey sphere). The star in the center of 
the sphere represents the MM state. 

state. The three decohered trajectories are different but equivalent, and a nice symmetrical 
configuration is observed. 

These results can be extended to robust states in higher dimensions and also to the GHZ 
and W states under the action of the BF and BPF channels. The results obtained by using 
the Hilbert-Schmidt distance instead of the QJSD are equivalent. The most representative 
distances that define the geometry shown in Fig. 3 are given in the following table: 



states 


djs 




initial-final 


0.6548 


0.9129 


initial-MM 


0.8285 


0.9682 


final-MM 


0.4188 


0.3227 


final-final 


0.6352 


0.4546 



As already mentioned, due to its highly symmetric character, the depolarizing channel 
does not have a single robust state. Any maximally entangled state will be robust for this 
decoherence process. All the associated trajectories in state space leading to these maximally 
entangled states are equivalent, and the sudden death of the entanglement always occurs for 
states with the same degree of mixedness. Systems undergoing this process evolve to a single 
asymptotic state, i.e., all initial states evolve asymptotically to the MM state. Moreover, 



10 



10 p i 1 i 1 s • • 1 • 1| 

0.8 - 

0.6 - ° o E - 

0.4 - 

o 

0.2 - 

O 

0.0 U ' 1 ' 2 o 6 o 6 o 6 

0.0 0.2 0.4 0.6 0.8 1.0 

P 

FIG. 5: Entanglement evolution (empty dots) and linear entropy of the time evolved state (filled 
dots) for 10"^ maximally entangled states obtained from \HS) by applying local unitary transfor- 
mations, under depolarizing channel. All plotted quantities are dimensionless. 

they do it following equivalent trajectories. Fig. [5] depicts the decay of the entanglement 
and the linear entropy of the decohered state for 10^ initial maximally entangled states 
(equivalent under unitary operations) due to the action of the depolarizing channel. 

IV. CONCLUSIONS 

In the present work we have investigated the state space trajectories associated to pre- 
viously determined four-qubit robust states undergoing different decoherence processes. We 
have characterized these trajectories by computing the distance between the time dependent 
state resultant from the decoherence process and some fixed states, namely the initial, final, 
and maximally mixed (MM) states. We have shown that, for states that are robust under 
decoherence, i.e., those with the maximal amount of entanglement during the evolution, 
the trajectories' aspect may look different. However, by reference to Fig. 4, one detects a 
significant degree of equivalence and symmetry among them. 

In the case of depolarizing decoherence processes the final state is the MM-one, but the 
system becomes separable before reaching it, i.e., a sudden death of entanglement takes 
place for this process. According to the symmetry of this channel, all maximally entangled 
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states are robust, and they evolve in equivalent fashion. Because of this equivalence, the en- 
tanglement's sudden death for any initial maximally entangled state is always characterized 
by the same degree of mixedness. 
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